Structure and physical properties of new layered iron oxychalcogenide BaFe20Se2 
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We have successfully synthesized a new layered iron oxychalcogenide BaFe2 0Se2 single crystal. 
This compound is built up of Ba and Fe-Se(O) layers alternatively stacked along the c-axis. The 
Fe-Se(O) layers contain double chains of edge-shared Fe-Se(O) tetrahedra that propagate along the 
b-axis and are bridged by oxygen along the a-axis. Physical property measurements indicate that 
BaFe2 0Se2 is a semiconductor without the Curie- Weiss behavior up to 350 K. There is a possible 
long range antiferromagnetic (AFM) transition at 240 K, corresponding to the peak in specific heat 
measurement and two glassy transitions at 115 K and 43 K. The magnetic entropy up to 300 K is 
much smaller than the expected value for Fe 2+ in tetrahedral crystal fields and Mosbauer spectrum 
indicates that long range magnetic order is unlikely at 294 K. Both results suggest that a short 
range magnetic correlations exist above the room temperature. 

PACS numbers: 75.50.Ee, 75.10.Pq, 75.30.Gw, 75.47.Np 



I. INTRODUCTION 

The discovery of layered iron oxypnictides 
(Ln(0,F)FeAs, Ln = rare earth elements, 1111-type) 
superconductors with T c up to 56 Kji has stimulated 
great interest in mixed- anion materials. Mixed anions 
from the same row of the Mendeleev periodic table 
tend to randomly occupy the same crystallographic site 
(anion disorder) because of the relatively similar sizes. 
On the other hand, if the anions are from different rows 
(oxysulfide or oxyselenide compounds for example), the 
distinctive difference of anion sizes and ionic polarization 
might lead to ordered occupancy in the different crys- 
tallographic sites (anion order), often forming layered 
crystal structure* 2 - 

The mixed-anion compounds have attracted some in- 
terest during the exploration of novel cuprate high 
temperature superconductors. Halooxocuprates are ex- 
ample where copper ions are coordinated with four 
oxygen ions, forming the Cu02 sheet, whereas the 
halogen ion usually occupies a so-called apical site. 3 
With proper electron or hole doping, halooxocuprates 
will become superconductors. 3 Besides superconductiv- 
ity, mixed-anion materials also exhibit diverse physical 
properties. Copper oxychalcogenides LnCuOCh (Ch = 
S, Se, and Te), isostructural to 1111-type iron-based su- 
perconductors, are wide-gap p-type semiconductors with 
transparent p-type conductivity, photoluminescence, and 
large third order optical nonlinearity^ - — Transition 
metal oxychalcogenides Ln20 2 TM 2 OCh2 (TM = Mn, Fe, 
and Co) contain the layers built up by the edge-shared 
octahedral unit [TM2OCI12] 2- . These materials show 
strong electron-electron interactions (Mott insulators) on 
the two dimensional (2D) frustrated antiferromagnetic 
(AFM) checkerboard spin-lattice^ - — 



The materials with [TM20Cli2] 2_ layers exhibit sim- 
ilar structural diversity to iron-based superconduc- 
tors. This is understandable since new compounds 
can be obtained by simply replacing [FeAs] - lay- 
ers with [TM2OCI12] 2 " layers, such as LnOFeAs -» 
Ln 2 2 TM 2 OCh2, AEFFeAs (AE = alkali earth metals) 
AE2F 2 TM 2 OCh2^"-2- 2 and AFeAs (A = alkali met- 
als) — >• Na2Fe20Se2^ _ — Thus, it is of considerable in- 
terest to explore the structural derivatives of AEFe2As2 
among the oxychalcogenide compounds. 

Here, we report the detailed synthesis and physi- 
cal properties of a new layered iron oxychalcogenide 
BaFe20Se2 single crystal. Even though it has the same 
chemical formula and [TM20Cli2] 2_ layers, the struc- 
ture is different from other iron oxychalcogenides with 
[TM 2 OCh 2 ] 2 " layers. To the best of our knowledge, 
BaFe20Se2 is the first layered iron oxychalcogenide with 
alkali earth metal. It shows a semiconducting behavior 
with possible successive spin-glass transitions at low tem- 
perature and short range antiferromagnetic order above 
the room temperature. 



II. EXPERIMENT 

Single crystals of BaFe20Se2 were synthesized by a 
self- flux method. Ba rod, Fe powder, Fe203 powder and 
Se shot were used as starting materials. BaSe was pre- 
reacted by reacting Ba piece with Se shot at 800 °C for 
12 hours. BaSe was mixed with other reagents and inti- 
mately ground together using an agate pestle and mortar. 
The ground powder was pressed into pellets, loaded in an 
alumina crucible and then sealed in quartz tubes with Ar 
under the pressure of 0.2 atmosphere. The quartz tubes 
were heated up to 600 °C in 10 h, kept at 600 °C for 12 
h, ramped again to 1100 °C in 12 h, kept at 1100 °C for 



12 h, and then cooled slowly at a rate of 3 °C/h to 800 
°C, finally the furnace was shut down and the sample 
was cooled down to room temperature naturally. Black 
plate-like crystals with typical size 2x2x0.5 mm 3 can be 
grown. Except for heat-treatment, all of processes for 
sample preparation were performed in glove boxes filled 
with argon. 

The crystal structure of BaFe20Se2 crystal was iden- 
tified by single crystal x-ray diffraction (XRD). The 
data were collected using the Bruker APEX2 software 
package^ 6 on a Bruker SMART APEX II single crystal 
x-ray diffractometer with graphite-monochromated Mo 
K a radiation (A = 0.71073 A) at room temperature. 
Data processing and an empirical absorption correction 
were also applied using APEX2 software package. The 
structures were solved (direct methods) and refined by 
full-matrix least-squares procedures on \F 2 \ using Bruker 
SHELXTL program software. 17 Obtained lattice param- 
eters of BaFe20Se2 are given in Table 1. Atomic coor- 
dinates, isotropic displacement parameters and selected 
interatomic bond distances and angles are listed in Table 
2. Phase identity and purity were confirmed by powder 
X-ray diffraction carried out on a Rigaku miniflex X-ray 
machine with Cu K a radiation (A = 1.5418 A). Structural 
refinements of powder BaFe20Se2 sample was carried out 
by using Rietica software. 18 

The average stoichiometry was determined by exami- 
nation of multiple points using an energy-dispersive X- 
ray spectroscopy (EDX) in a JEOL JSM-6500 scanning 
electron microscope. The presence of oxygen was con- 
firmed for BaFe20Se2, but the exact amount could not 
be quantified because of experimental limitations. The 
average atomic ratios determined from EDX are Ba:Fe:Se 
= 1.00(7):2.2(2):2.0(2), close to the ratio of stoichiomet- 
ric BaFe20Se2- 

Mosbauer spectrum was taken in transmission mode 
with 57 Co(Rh) source at 294 K and the parameters were 
obtained using WinNormos software. 20 Calibration of the 
spectrum was performed by laser and isomer shifts were 
given with respect to a-Fe. 

The X-ray absorption spectra of the Fe and Se if -edges 
were taken in transmission mode on powder samples of 
BaFe20Se3 at the X19A beamline of the National Syn- 
chrotron Light Source. Standard procedure was used to 
extract the extended x-ray absorption fine-structure (EX- 
AFS) from the absorption spectrum^ 

Electrical transport, heat capacity, and magnetiza- 
tion measurements were carried out in Quantum Design 
PPMS-9 and MPMS-XL5. 



III. RESULTS AND DISCUSSION 

Similar to K x Fe2- 2/ Se2 and BaFe2Se3 j 21 i 22 the struc- 
ture of BaFe20Se2 is built up by stacking the Ba cations 
and Fe-Se(O) layers alternatively along the c axis (Fig. 
1(a)). However, within the Fe-Se(O) layers, the connec- 
tion of Fe-Se(O) tetrahedra in BaFe20Se2 is distinctive 
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FIG. 1. (a) Crystal structure of BaFe2 0Se2. The biggest 
orange, big red, medium blue and small purple balls represent 
Ba, Fe, O and Se ions, respectively, (b) and (c) Powder and 
single crystal XRD patterns of BaFe2 0Se2. 



different from those in K x Fe2- 2/ Se2 and BaFe2Se3. In 
BaFe20Se2, double chains of edge-shared Fe-Se(O) tetra- 
hedra propagate along the Se atoms parallel to b- axis. 
The Fe-Se(O) double chains are bridged by oxygen along 
the a-axis (Fig. 1(a)). The bond distances of Fe-Se (~ 
2.5 A) andFe-0 (~ 1.9 A) are very similar to distances in 
compounds where Fe is also in tetrahedral coordination 
(BaFe 2 Se 3 : d Fe - S e ~ 2.4 A, Fe 3 4 : d Fe -o ~ 1.9 A). Be- 
cause of different bond distances between Fe-Se and Fe-O, 
Fe atoms are located at the highly distorted tetrahedral 
environment when compared to other materials with pure 
Fe-Se tetrahedra. This distortion is also reflected in the 
significant deviation of bond angles from the value for 
the ideal tetrahedron (109.5°). The bond angles range 
up to 130.38(10)° in BaFe20Se2, much larger than the 
values in BaFe2Se3 and K x Fe2- 2/ Se2^^ Moreover, the 
different connection of Fe-Se (O) leads to the larger near- 
est neighbor Fe-Fe distances (dFe-Fe (~ 3.13 A)) when 
compared to BaFe2Se3 and K iC Fe2- 2/ Se2^^ Using ob- 
tained Fe-Se and Fe-0 bond lengths as shown in Table 
2, the valence of Fe ions in BaFe20Se2 can be calcu- 
lated using the bond valence sum (BVS) formalism in 
which each bond with a distance dij contributes a va- 
lence Vij = exp[(i?ij — dij)/ 0.37] with R^ as an empir- 
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TABLE I. Crystallographic Data for BaFe 2 OSe 2 . 



Chemical Formula 
Formula Mass (g/mol) 
Crystal System 
Space Group 
a (A) 
b(A) 

c(A) 
V (A 3 ) 
Z 

Density (g/cm 3 ) 

Rl/wR2 



BaFe 2 OSe 2 
422.94 
orthorhombic 
Pmmn (No. 59) 
9.8518(7) 
4.1332(3) 
6.7188(4) 
273.59(3) 
2 

5.134 



F > 4aF )/Rl (all data) a 0.0459/0.1423/0.0476 
Goodness-of-Fit 1.187 

a R1 = Z\F \-\F C \/Z\F \, 
wR2 = p(|F 2 - I^DVS^Fq 2 ) 2 ] 1 / 2 . 
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velocity (mms 1 ) 

FIG. 2. Mosbauer spectrum of BaFe 2 OSe 2 at T = 294 K. Inset 
shows local coordination of Fe atoms in tetrahedra: central 
smallest ball denote Fe, medium size ball shows oxygen wheras 
largest balls show selenium atoms. 



ical parameter and the total of valences of atom z, Vi 
equals Vi = ^2 v ij The calculated valence of Fe ions 

3 

is +2.23, slightly larger than the apparent oxidation state 
(+2) for Fe ions. Structural refinement of powder XRD 
results confirms that all reflections can be indexed in the 
Pmmn space group and the refined lattice parameters 
are a = 9.862(1) A, b = 4.138(1) A, and c = 6.732(1) 
A, close to the values obtained from single crystal XRD. 
Only (001) reflections were observed in the single crys- 
tal XRD patterns of a BaFe20Se2 (Fig. 1(c)), indicating 
that the crystallographic c axis is perpendicular to the 
plane of the single crystal. 

Mossbauer fit of BaFe20Se2 consists of two doublets 
corresponding to the paramagnetic state (Fig. 2). Based 
on the goodness of fit (0.9995) long range magnetic or- 
der is unlikely at 294 K. The first doublet dominates the 
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FIG. 3. FT magnitudes of the EXAFS oscillations (symbols) 
for Fe K-edge (a) and Se K-edge (b). The model fits are shown 
as solid lines. The FTs are not corrected for the phase shifts 
and represent raw experimental data. Insets of (a) and (b) 
filtered EXAFS (symbols) with k-space model fits (solid line). 



signal and takes 92 % of the fit area with isomer shift 
of 0.690(1) mm/s. First coordination sphere of Fe (in- 
set of Fig. 2) is tetrahedral with three selenium and one 
oxygen atoms. Fe 2+ and Fe 3+ in Fe04 tetrahedra have 
isomer shift values (0.85-0.95) mm/s and (0.2-0.3) mm/s 
respectively. Isomer shift drops with the increase of ef- 
fective charge near 57 Fe. Smaller value of isomer shift in 
BaFe20Se2 when compared to Fe04 tetrahedra is consis- 
tent with larger covalency of Fe-Se bonds. Isomer shift 
is in agreement with the calculate value for Fe valence 
by BVS formalism. Quadrupole shift of 1.594(2) mm/s 
points to large distortion of tetrahedra. Second doublet 
has isomer shift of 0.51(1) mm/s and 0.76(3) quadrupole 
splitting. These parameters are similar to what was ob- 
served in other Fe compounds (FeSe2 or a-FeSe), point- 
ing either to presence of impurities or to intrinsic prop- 
erty of Fe environment in BaFe20Se2- Assuming similar 
recoil-free factors 8 % of fit area corresponds to approxi- 
mately similar volume (crystallographic) phase fraction. 
Since such large impurity fraction would have been de- 
tected in our X-ray experiments this might suggest that 
second doublet comes from undistorted FeSe4 tetrahedra 
in BaFe20Se2 with no oxygen bonds present. In other 
words, there might be 8% ionic disorder existing. 

The local environments of Fe and Se such as bond 
distances are also analyzed by fitting the EXAFS data 
as show in Fig. 3. For the Fe site, the nearest neigh- 
bors are one oxygen atom (d Fe - = 1.8769(19) A) and 
three Se atoms (lxSe(I) and 2xSe(II)) with two differ- 
ent distances (d Fe _se(i) — 2.4706(6) A and d Fe _s e (ii) — 
2.5540(4) A). The next nearest neighbors are three Fe 
atoms (lxFe(I) and 2xFe(II)) with two different dis- 
tances (^Fe-Fe(I) ~ 3.13 A and d F e-Fe{II) ~ 3.24 A). 

On the other hand, for the Se site, the nearest neigh- 
bors are three Fe atoms with bond distances d Fe _s e (i) 
and d Fe _se(ii)- The next nearest neighbors are three 
Ba atoms (2xBa(I) and lxBa(II)) with two different 
distances (d Se -Ba(i) = 3.3541(3) A and d Se -Ba{ii) = 
3.4008(4) A). From the joint analysis of Fe and Se edges 
EXAFS data using a single bond distance for Fe-O, Fe- 
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TABLE II. Atomic Coordinates, Equivalent Isotropic Displacement Parameters, and Selected Bond Lengths and Angles for 
BaFe 2 OSe 2 . 



Atom 


W y CKOti 




y 


z u eq (A 2 ) 


oa 


ZD 


3/4 


1/4 


0.49047(4) 0.00734(9) 


re 


A £ 
41 


0.58580(5) 


3/4 


0.87890(7) 0.00915(11) 





2a 


3/4 


3/4 


0.7372(6) 0.0089(5) 


Se 


4f 


0.45734(4) 


1/4 


0.75887(5) 0.00739(10) 






Interatomic Distances (A) 


Fe-Se 


2.4706(6) 


Fe-Se 


2.5540(4) 


Fe-0 1.8769(19) 


Ba-Se 


3.3541(3) 


Ba-Se 


3.4008(4) 


Ba-0 2.649(2) 


Ba-Fe 


3.7012(4) 









Bond Angles (°) 



Se-Fe-Se 103.050(15) Se-Fe-Se 108.03(2) Se-Fe-0 130.38(10) 
Se-Fe-0 105.48(5) Fe-Se-Fe 76.950(15) Fe-Se-Fe 108.03(2) 
Fe-O-Fe 119.1(2) 
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FIG. 4. Temperature dependence of the resistivity p a b(T) of 
the BaFe2 0Se2 crystal with H = (closed blue square) and 
90 kOe (open red circle, H||c). Inset shows the fitting result 
of pabiT) at zero field using thermal activation model. The 
red line is the fitted curve. 



Se, Fe-Fe and Se-Ba, and by fitting the k range 2-14 
A -1 for Fe if -edge and 2-12.9 A -1 for Se if -edge (insets 
of Fig. 3(a) and (b)), the fitted average bond lengths 
are d Fe -o = 1.87(2) A, d Fe - S e = 2.500(8) A, d Fe - Fe = 
3.17(4) A and dse-Ba = 3.38(2) A, which are consistent 
with the average bond distances derived from XRD fit- 
ting. 

As shown in Fig. 4, the resistivity p a b(T) of BaFe20Se2 
crystal increases rapidly with decreasing temperature, 
suggesting that this compound is a semiconductor in 
measured temperature region. The semiconducting be- 
havior of pab(T) can be fitted using the thermal acti- 
vation model PabiT) = po exp(i? a //c£T), where po is a 
prefactor, E a is thermal activated energy and ks is the 
Boltzmann's constant (inset in Fig. 4). The obtained 



E a is 0.360(2) eV and the room-temperature resistivity 
p a b(300K ) is about 10 k^-cm. Both are much larger than 
in BaFe 2 Se 3 (0.178 eV and 17 ^-cm). 22 The E a is also 
larger than in La202Fe20Se2 where Fe has octahedral 
coordination with 4xSe and 2x0 mixed anions *2£ The 
large p a b(300K) and E a could be due to the shorter Fe- 
O bond in highly distorted structure which localizes elec- 
trons and increases the band gap. Since La202Fe20Se2 
is a Mott insulator with narrow 3d electronic bands due 
to strong correlation effects, it is of interest to study 
whether BaFe20Se2 is normal band insulator or Mott 
insulator via theoretical calculations. The p a b(T) mea- 
sured at H — 90 kOe indicates that there is no obvious 
magnetoresistance in BaFe20Se2 (Fig. 4). 

The Curie- Weiss temperature dependence of x(T) is 
not observed in BaFe20Se2 single crystals up to 350 K 
(Fig. 5(a)). On the other hand, Mossbauer fit indicate 
that long range magnetic order is absent at 294 K. It 
implies that short range magnetic interactions might be 
present already at the room temperature. There is a peak 
in the d(xT)/dT curve at T x = 240 K, 27 for H|| and the 
values of susceptibility between ZFC and FC curves are 
nearly the same at this temperature. This peak corre- 
sponds to the AFM transition temperature T/v = 240 K, 
confirmed by specific heat measurements. Moreover, the 
decrease in Xc(T) with temperature is more significant 
than Xab(T). It suggests that the easy- axis of magneti- 
zation direction could be in afr-plane. According to mean- 
field theory for collinear antiferromagnets, x(T) — >• 
along the easy-axis direction whereas it is nearly constant 
below T/v for the field perpendicular to the easy-axis di- 
rection. With further decreasing temperature, there are 
two peaks in the d(xT)/dT curve at T2 = 115 K and T3 
= 43 K, respectively (inset of Fig. 5(a)). The hysteresis 
between zero-field-cooling (ZFC) and field-cooling (FC) 
measurements implies that they are spin-glass-like tran- 
sitions. The absence of hysteresis in isothermal M(H) 
loops for H || ab (Fig. 4(b)) at T = 1.8 K and 250 
K indicates that there are no ferromagnetic impurities. 
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FIG. 5. (a) Temperature dependence of DC magnetic suscep- 
tibility x(T) under ZFC and FC modes with H = 1 kOe along 
ab-plane and c-axis. (b) Isothermal magnetization hysteresis 
loops M(H) for H\\ab at T = 1.8 K and 250 K. 



Moreover, the slope of M(H) increases with increasing 
temperature consistent with the AFM behavior observed 
in x(T) curves (Fig. 5(b)). 

As shown in Fig. 6, the specific heat C p of BaFe20Se2 
approaches the value of 3NR at high temperature (T = 
300 K) assuming N = 6, where N is the atomic number 
in the chemical formula and R is the gas constant (R = 
8.314 J mol -1 K _1 ). It confirms the atomic numbers of 
BaFe20Se2 obtained from single crystal XRD fitting and 
consistent with the atomic ratio measured from EDX. On 
the other hand, at the low temperature, C P (T) curve can 
be fitted by a cubic term f3T 3 (inset of Fig. 6) and the 
fitted value f3 is 0.852(2) mJ mol -1 K -4 . According to 
the formula 9^ = (127r 4 iVi?/5/?) 1/3 , Debye temperature 
is estimated to be 6^ = 239.1(2) K. The larger when 
compared to BaFe2Se3 (Q D = 205(1) K) can be ascribed 
to the smaller atomic mass of oxygen than selenium.— 

There is a A- type anomaly at Ti = 241.7 K as shown in 
Fig. 6. The peak position is very close to the peak posi- 
tion (Ti = 240 K) at dx(T)/dT curve (inset of Fig. 5(a)). 
It suggests that a long-range AFM ordering forms at this 
temperature. After substraction the phonon contribution 
(C p h) fitted using a polynomial for the total specific heat, 
the magnetic contribution (C ma g) can be obtained. The 
magnetic entropy can be calculated using the integral 




FIG. 6. 



150 
T(K) 

Temperature dependence of specific heat for 
BaFe2 0Se2 crystal. The horizontal orange line represents the 
classical value according to Dulong-Petit law at high tem- 
perature with N = 6. The blue curve represents the phonon 
contribution fitted by a polynomial. The right axis and its as- 
sociated green solid curve denote the magnetic entropy. Inset: 
the low- temperature specific-heat data in the plot of C p /T vs 
T 2 . The red solid line is the fitting curves using formula 
Cp/T f3T 2 . 



Sma 9 (T) = ^ C mag /Tdt. The derived S ma? is ~ 1.83 
J/mol-K when T is up to 300 K (Fig. 6), which is much 
smaller than expected values Rln(2S' + 1) for Fe 2+ ions 
with low spin state (~ 20% Rln3) and with high spin state 
(~ 13.7% Rln5) in tetrahedral crystal fields. It is even 
less than Rln2 for S=l/2 (~ 31.8% Rln2), suggesting that 
there is short range magnetic order at higher tempera- 
ture which partially releases the magnetic entropy before 
any long range magnetic transition, consistent with the 
Mossbauer spectrum and magnetization measurement re- 
sults. On the other hand, there is no peaks appear at T3 
= 115 K and T4 = 43 K, confirming the glassy magnetic 
behaviors at those temperatures. 



IV. CONCLUSION 

In summary, we report a discovery of a new layered 
iron oxychalcogenide BaFe20Se2- The crystal structure 
features double chains of edge-shared Fe-Se(O) tetra- 
hedra that propagate along the 6-axis. The Fe-Se(O) 
double chains are bridged by oxygen along the a- axis. 
BaFe20Se2 single crystals are magnetic semiconductors 
with a possible long range AFM transition at 240 K and 
two glassy transitions at 115 K and 43 K. The magnetic 
entropy up to 300 K is much smaller than expected for 
Fe 2+ in tetrahedral crystal fields, suggesting the existence 
of short range AFM correlation above room tempera- 
ture. Because of the interesting structure and connec- 
tivity among the iron atoms, it is of interest to investi- 
gate the mechanism of magnetic ground state by theoret- 
ical calculations and by neutron scattering experiments. 
Moreover, doping effects on physical properties should 
further reveal similarities and differences with supercon- 
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ducting iron selenide compounds. 
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